2, 2008 1:30 WSPC/INSTRUCTION FILE VitevDPF2004 



International Journal of Modern Physics A 
© World Scientific Publishing Company 



PROBING THE PHASES OF QCD 
IN ULTRA-RELATIVISTIC NUCLEAR COLLISIONS 

IVAN VITEV 

Los Alamos National Laboratory, Theory Division and Physics Division 
Mail Stop H846, Los Alamos, NM 87545, USA 

Received 25 November 2004 

The status of RHIC theory and phenomenology is reviewed with an emphasis on 
the indications for the creation of a new deconfined state of matter. The critical role of 
high energy nuclear physics in the development of theoretical tools that address various 
aspects of the QCD many body dynamics is highlighted. The perspectives for studying 
nuclear matter under even more extreme conditions at the LHC and the overlap with 
high energy physics is discussed. 
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1. Introduction 

Quantum Chromodynamics and asymptotic freedonP predict the existence of a new 
state of matteiEI, the quark-gluon plasma (QGP), at exceedingly high temperatures 
and energy densities similar to the ones that characterized the first few microseconds 
after the Big Bang. The quest for such a deconfined state of QCD has led to the 
highly successful program at the Relativistic Heavy Ion Collider (RHIC) and will 
be an integral part of the program at the Large Hadron Collider (LHC). Recent 
theoretical developments have improved our understanding of the complex many- 
body dynamics in high energy nuclear collisions. Comparison between data and 
theory is suggestive of the creation of a deconfined state of QCD with energy density 
on the order of 100 times normal nuclear matter density. For a review of the RHIC 
experimental results see'^. 

2. Soft and Intermediate px Hadrons 

An economical description of the bulk particle production in A + A reactions can 
be achieved in the framework of the thermal model. All measured particle ratios 
at midrapidity from the y/sJ^N — 130, 200 GeV Au + Au runs at RHIC are well 
reproduced with Tj k.Tc~ 175 MeV and = 40, 30 MeV, respectiveljEI. Such 
approach, while instructive, does not carry information about the dynamical evo- 
lution prior to freeze-out. Constraining the initial conditions, relevant at the early 
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stages of relativistic heavy ion coUisions, requires the development and application 
of microscopic models. 

Relativistic hydrodynamics. Differential particle distributions at low px can be 
calculated in the framework of relativistic hydrodynamicJ^, which assumes local 
thermal equilibrium and solves the energy-momentum and current conservation 

d^T^^x)^0, d,j^ix)^0. (1) 

Additional constraints, necessary to determine the system Eq. include, for 
example, ideal hydrodynamics and equation of state p ~ p(e) or a boost invari- 
ant Bjorken model. Good description of the pt < 2 — 3 GeV hadron spectra for 
TT^, K^, p{p), A(A), , • • • can be obtained. However, the elliptic flow V2 of mas- 
sive hadrons was shown to be much more sensitive to the choice of equation of state 
(EOS) with first order phase transition (preferred) versus hadronic EO^. Traces 
of early partonic thermalization, as suggested by hydrodynamics, can also be found 
in the in the mean pT fluctuations versus centralitj. . 

Recent hydrodynamic simulations have employed initial conditions motivated 
by gluon saturation phenomenology. . In A + A reaction this approach alleviates 
the problem with the transverse energy and the parton number in such models. 
However, it leaves an open question for p + A reactions where the hydrodynamic 
description is not applicable. 

Parton Coalescence. Another microscopic approach to low and moderate pt 
hadroproduction is the covariant Boltzmann transporl^. T he t heory was shown to 
be applicable not only in A + A but also in p + A reaction^^ and has provided a 
good description of the rapidity density of hadrons dN^^ /dy and their transverse 
momentum distributions. The output of a partonic transport can be incorporated in 
the coalescence/recombination models that convolute the quark Wigner functions 
with the meson and baryon wavefunctions"'^"'^. These predict scaling properties of 
the elliptic flow v 

v^\pt) = 2vI{pt/2) , vUpt) = HiPT/i) (2) 

at intermediate transverse momenta and enhanced baryon to meson ratios such as 
p{p)/n^ and A{A)/K^. The energy loss of quarks and gluons plays an important 
role in the generation of the elliptic flow at the partonic levejJ>2« It also suppresses 
the perturbative high pr hadron production and thus extends the region of possible 
non-perturbative contribution for baryon^^l. Coalescence models have been used 
to calculate the expected charm meso n tr ansverse momentum distributions and the 
charm meson and baryon elliptic flo'vA^. 



3. The Perturbative QCD Factorization Approach 

In the nuclear environment, modifications to the QCD factorization approacfU^ 
arise from the elastic, inelastic and coherent multiple scatte ring . These can be sys- 
tematically incorporated in to the perturbative formalisniLt' For £ + A, p + A 
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Fig. 1. Left panel frorJi^ resummed QCDpower corrections to nuclear shadowing in £ + A, 
manifest in F2{A)/ F2{D). Right panel frorrl. modification of the away side correlation function 
C(A(/)) versus rapidity, centrality and px from higher twist. 



and A + A reactions their effect on experimental observables is summarized in 
Table 1. The relative importance of these distinct contributions in the different 
kinematic and density regimes is yet to be determined. Elastic scattering leads to a 
Cronin enhancement and a small broadening of the away side cor relation function 
C(A0) = {l/Ntrig)dN''^''^/dA(t) . These have been discussed irPl. 

Coherent power corrections. A class of corrections that can be naturally incorpo- 
rated in the perturbative QCD factorization approach is associated with the power 
suppressed ~ n > 1 contributions. The higher twist terms are typically ne- 

glected in reactions with "elementary" nucleons for > 1 GeV^. However, in the 
presence of nuclear matter such corrections can be enhanced by the large nuclear 
size ^ A^/^. 

The effective longitudinal interaction length probed by the virtual meson of 

Table 1. Effect of clastic, inelastic and coherent multiple scattering on the transverse momentum dis- 
tribution of hadrons in the perturbative regime. 

Type of scattering Transverse momentum dependence of the nuclear effect 



Elastic 

Inelastic 

Coherent 



Suppression at low px, enhancement at moderate pT, disappears at high pT 
Suppression at all px, weak px dependence, persists at high px 
Suppression at low pT, disappears at high pT 
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momentum (7^ is characterized by l/xP. If the momentum fraction of an active 
initial-state parton x ^ Xc — l/2mArro ~ 0.1 with nucleon mass m^r and ra- 
dius tq, it could cover several Lorentz contracted nucleons of longitudinal size 
~ 2ro{rnN / P) in a large nucleus. Each of the soft interactions is characterized by a 
scale of power correction per nucleon « '^^g'l'^'* {p\F'^\p) with the matrix element 

{p\F^\p) = ^lYmx^oxG{x,Q^]^^. The multiple final state scattering of the struck 
quark generates a dynamical parton mass ^^yn ~ — 1) and a consequent 

rescaling in the value of Bjorken-a^lSl 

XB XB + XB-^ + XB ^2 =a;slH ^2 I • 

In Eq. Q is the physical mass of the quark in the final state. 

In Ref. 1^ next-to-leading order global analysis of the nuclear parton distribu- 
tions (nPDFs) was shown to strongly disfavor more than 5% gluon shadowing. This 
result also suggests that the observed attenuation in the structure functions is a 
result of multiple final state scattering since any dynamical mechanism will predict 
an effect twice as large for gluons when compared to quarks, — > Ca/Cf^^- 

Application to nuclear shadowing in the DIS structure function F2 {A) is shown in 
the left panel of Fig. 1. The coherent nuclear enhanced power corrections in single 
a nd double inclusive hadron production m p + A reactions have been evaluated 
The apparent suppression of the away side correlation function C{/S.(f)) = 
(1 / Ntrig)dN^^^'^ / dlS.(f) is shown in the right hand side of Figure 1. As seen from Eq. 
0, these modifications disappear at high pT- It will be instructive to investigate 
the open charm production^ including resummed power corrections. Additional 
discussion of n uclear shadowing and antishadowing in the light front approach can 
be found irPS. 

Non-Ahelian energy loss. The most efficient mechanism that modifies the large 
transverse momentum hadroproduction in the presence of a hot and dense quark- 
gluon plasma is the medium induced gluon bremsstrahlung . It leads to strong 
suppression of high-p^ particles, known as jet quenching. Non-Abelian energy loss 
was also shown to be an effective attenuation mechanism in cold nuclear matteJ^. 

The qualitative behavior of the energy loss as a function of the density and the 
size of the system can be calculated using the GLV approacfP^. To first order in 
opacity for static and 1-l-lD Bjorken expanding plasmas 

(A£> J"^"'™''"^- (4) 

ApplicatiorPS of the QCD theory of energy loss to the single inclusive pion sup- 
pression and the attenuation of the away side dijet correlations is shown in Fig. 2. 
These persist, as shown by the calculation, to high pT- For similar results see^. 

Better experimental techniques may reveal the redistribution of the lost jet 
energj^' in low frequency modes > ujpi. The quenching of charm and beauty 
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Fig. 2. Left panel fron® nuclear modification factor Ra_a{pt) S't ^sjvjv = 17, 62 and 200 GeV. 
Right panel fronJ^l modification of the away side correlation function C(A(/)) from energy loss at 
high Pt- 



mesonJ^will provide a physical scale, and mt, relative to which the parameters 
of the medium can be reliably extracted. 



4. Conclusions 

The complex dynamics of ultrarelativistic heavy ion reactions invites a variety of 
theoretical approaches, applicable in the different kinematic and density regimes, in 
the description of experimental observables. Consistent application of the models is 
required to better understand their strengths and limitation and to disentangle the 
distinct nuclear effects. At present, combined results from relativistic hydrodynam- 
ics and jet tomography suggest the creation of QGP at RHIC with initial energy 
density e ~ 15 — 20 GeV/fm'^. The interactions in the plasma, however, remain 
strong and are consistent with the lattice expectation of as — 0. m 

The upcoming RHIC runs and the LHC program will ensure the dominance 
of high Pt and jet physicJ^ and will provide precision A + A data at ^ysNN = 
200 - 5500 GevEHl. This opens new possibilities for studying photon and dilepton 
tagged jets, heavy flavor and multiparticle correlations. For cold nuclear matter, a 
low energy p{d) + A run at RHIC will be a critical step in clarifying the relative 
importance of the nuclear effects summarized in Table 1. 
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